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1 Introduction 

In the last decade chaotic behaviour in general relativity started to be the 
subject of many interesting papers. Two main lines of research can be recog- 
nized. The first deals with chaoticity associated with inhomogeneous cosmo- 
logical models, the second line assumes a given metric and looks for chaotic 
behaviour of geodesic motion in this background. An interesting selection of 
references can be found in . In particular many papers devoted to the study 
of chaotic dynamics in general relativity resort to the Poincare- Melnikov 
method (see e.g. [0, |3| for the general theory). The Melnikov method is 
an analytical criterion to determine the occurrence of chaos in integrable 
systems in which homoclinic (or heteroclinic) manifolds biasymptotic to un- 
stable critical points or to periodic orbits (more generally to invariant tori) 
are subjected to small perturbations. Such perturbations may lead to the 
phenomenon of transversal intersections of the stable and unstable manifolds. 
This kind of dynamics can be then detected by the Melnikov functions, since 
they describe the transversal distance between the stable and the unstable 
manifolds of the critical point or periodic orbit. 

The Melnikov method has been applied in many branches of physics and 
applied mathematics, so that it is impossible to give here even a partial 
account of the vast literature, and we will quote only some of the applications 
to general relativity that are more strictly related to the present paper (for a 
more complete list of references see e.g. |], ]5[]). Examples of applications 
of the Melnikov method in general relativity concern the study of the orbits 
around a black hole perturbed either by an external quadrupolar shell 0, || 
or gravitational radiation [1], [7| . 

In this work we take a slightly different approach, compared with the 
previous literature. Firstly, we do not analyze perturbations of the metric, 
but perturbations due to the interactions produced by uniform electric or 
magnetic fields. We consider these as perturbations to the Hamiltonian of 
a charged test particle in free fall in a Schwarzschild black hole. Secondly, 
we do not restrict ourselves to perturbations lying in the plane of the orbit, 
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but we deal more in general with the full three dimensional problem, i.e. 
with perturbations which may change the plane of the orbit. Let us also 
remark that we will consider time-independent perturbations, and show that 
they produce chaos; this happens even when only the rotational invariance 
is broken. 

The problem we study in this paper can be related to the astrophysical re- 
ality, since it is well known that magnetic fields can be associated with black 
holes. We will consider Schwarzschild black holes, i.e. non rotating ones, 
and we reserve to tackle the more interesting, but more complex problem of 
rotating black holes in a forthcoming paper. Isolated black holes cannot pos- 
sess any properties other than mass, electric charge and angular momentum, 
but the medium surrounding the hole can be responsible for the magnetic 
field. For example, super massive black holes can acquire surrounding mat- 
ter either by gravitationally pulling interstellar gas into its vicinity, or by 
the disruption of passing stars. The surrounded matter will be shaped in an 
accretion disk in a state of plasma, that produces a magnetic field. 

Electric fields are less likely to be found near a black hole, although the 
Blandford-Znajek mechanism is a process that can develop a potential 
difference (i.e. an electric field) between the poles and the equator of a black 
hole spinning in a magnetic field pointing along the axis of rotation. 

On the other hand, it is not easy to account for the presence of charged 
particles since the ionized gas, that is shaped in an accretion disk surrounding 
a black hole, is in a state of plasma, and a plasma is electrically neutral. One 
way to account for the presence of charged particles is to consider again the 
Blandford-Znajek mechanism, that, if the field strength is large enough, can 
separate charges and accelerate them to relativistic velocities. Otherwise, 
we can model the motion of electromagnetic currents in a macroscopic piece 
of plasma considering the equivalent problem of the motion of a charged 
particle moving around a black hole surrounded by a magnetic field. This 
can be done since the magnetic field has the same effect on an electric current 
in a macroscopic piece of plasma as on a single charge. 

As said before, we will deal in this paper only with uniform fields. Al- 
though the fields which can be found in the astrophysical reality, e.g. the 
Blanford-Znajek ones, are not uniform, some authors (see || and references 
therein) pointed out that models with uniform external fields are a fairly 
good approximation in order to explain the qualitative features of a black 
hole in a magnetic field, indeed qualitative arguments indicate the existence 
of a quasi-uniform field. 
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In the next section we discuss the equation of motion for the Schwarzschild 
solution, we find the homoclinic orbit and show that the perturbed systems 
we are considering are of type III, according to the classification given in . 
In section 3 we present a summary of the Melnikov method for a system of 
this type. In section 4 we consider the perturbations given by uniform electric 
and magnetic fields. In the last section we apply the Melnikov method and 
we prove the occurrence of chaos in the perturbed system both in the case 
of electric and of magnetic field, but with a great difference: the occurrence 
of chaos is a first-order effect (in the strength of the field) in the case of the 
electric field, whereas it is a second-order effect in the case of the magnetic 
field. In particular, the Hamiltonian perturbed by a constant magnetic field 
turns out to be integrable at the first-order, and the first non-vanishing con- 
tribution to the Melnikov integral comes from the second-order term. On the 
other hand, the occurrence of chaos in this situation has been already proven 
T0|| by means of numerical arguments based on the study of the trajectories 



in the Ernst metric; it can be interesting to point out that our approach con- 
firms this result using a completely different (analytic) approach. As another 
interesting result, we obtain that in both cases (electric and magnetic fields) 
only the components of the fields on the plane of motion are responsible for 
the chaotic behavior, whereas the components normal to the plane are not. 

It can be significant to remark that, although the occurrence of chaos 
can often be expected since the equations of motion are integrable only in 
very special cases, the consequences of chaos may be greatly relevant and 
sometimes far-reaching. Indeed, we can observe that chaotic systems may 
exhibit a very rich dynamics, including also regions of stability, periodic 
orbits, regions of ergodicity and so on, that can also lead to macroscopically 
observable phenomena. We recall, just to mention few examples, that chaotic 
defocusing of light might make black holes bright [^T|] , (and hence observable) 
and that the depletion of the outer asteroidal belt in the Solar System might 
be explained with the existence of a chaotic region |T2fl . Moreover, this last 



example also reveals the richness of chaotic dynamic since it was shown that 
the only few asteroids found in the outer belt are resonant and dynamically 
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2 Equations of Motion 

In this section we study the motion of a relativistic (charged) particle in free 
fall in a gravitational field. To describe the motion of a particle moving in 
a space-time with a metric g a b it is convenient to consider the action (with 
c=l): 

m t • „ m f b 



S=- J ds or S = - J g ab x a x b ds (1) 

where x a (s) denotes the world-line of a particle and m denotes the mass. 
Using the same notations as in ffl], the canonical conjugate momentum to x a 
is p a = mg ao x b and satisfies the mass-shell constraint g ab p a Pb = —m 2 . The 
Hamiltonian of the system can be defined as: 

H = (2) 

Let our background metric be the Schwarzschild metric, i.e. the metric of a 
non-rotating black hole: 

ds 2 = -fdt 2 + f' l dr 2 + r 2 {d6 2 + sin 2 9d<j) 2 ), (3) 

where / = 1 — 2M/r. Then the Hamiltonian (2) becomes: 

„ i (^ + 4 + ^U-r'(-y (4) 

2m \ r l r l sin J \m J 



where 



and 



x dr 2 d9 
Pr = fnf — , p e =mr — (5) 
ds ds 



E = —p t = mf— , pj, = mr 2 sin 2 9^- and L 2 = p 2 e -\ — (6) 
ds ds sin 9 

are conserved quantities. We can introduce the following Hamiltonian, more 
suitable for our purposes (i.e. for finding the orbits of a particle in free fall 
in the Schwarzschild metric): 

E 2 -m 2 1 / 2 p 2 e p\ \ Mm Mp 2 Mp 2 

H ° = O = n~ \ P r + ~2 + 9 ■ 2a 3 3 ■ 2 a ' ^ ' 

2m 2m \ r z r A sin 9 1 r mr 6 mr 6 sin 9 
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which can also be written in terms of H.o as: 

2m r 

This new Hamiltonian defines a problem similar to the one discussed in [0], 
but in our case we study the motion of a particle in the three dimensional 
space. Exactly as in jl], we can show that our Hamiltonian admits an un- 
stable circular orbit 7 in the plane 9 = tt/2, together with an homoclinic 
loop (or, to be more precise, a one parameter family of homoclinic loops) 
biasymptotic to this orbit; the spherical symmetry of the problem implies 
that the same happens for every plane for the origin. 

Choosing for simplicity the plane 9 = 71/ '2, it can be easily shown that 
this unstable circular orbit 7 has radius 

(9) 



l + (3' 



where j3 = y 1 — 12M 2 m 2 /p 2 b < 1, and that there are homoclinic orbits to 

this invariant set for < (3 < 1/2. The equations of motion for the homoclinic 
orbits are: 



dr _^./ 2 (E 2 -m 2 _V\ Mm Mp|\ d4 = vl_ nm 



ds Urn \ 2m r 2 r mr 3 J ds mr 2 ' 

where the sign — (resp. +) holds for s < (resp. s > 0). Let us denote by 

R = R(s) and $ = $(s) (11) 

the expressions of r and <fi with respect to s for the homoclinic orbit with 
R(0) = r max (the turning point) and $(0) = 7T (i.e. the orbit having its 
axis coinciding with the and the point r u in the positive x). We 

can observe that the maximum value of r, along this homoclinic solution, is 
r max — 6M/(1 — 2/3), but we will not need the explicit expression for those 
functions, that can be found in (1[| (at least for R = R(s)), we only retain, 
as in [14, 15], the information that R(s) is an even function and $(s) is an 
odd function of s. 

In the following sections, we will study time-independent perturbations 
that destroy the spherical symmetry of the unperturbed system. This sit- 
uation is another example, in a completely different setting, of symmetry- 



breaking perturbations that were already examined in [14, 15] 
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The perturbations that we consider in this work can be written in Hamil- 
tonian form as 

H = H o + eW(r,<p,0,p r ,Pi,pe) (12) 

We want to show that the Hamiltonian in (12) defines a system of type III 
(according to the classification given in |J) with n = 1 and m = 2, so that we 
can apply the Melnikov theory developed for this kind of problems to prove 
that chaos occurs in the perturbed problem. 

According to ||, systems of type III can be written in the following 
general form, where we have denoted by ui (to avoid any confusion with the 
azimuthal variable 0) the angle (cyclic) variables conjugated to the action 
variables J: 

x = JD x H (x, I) + eJD x W(x, I, u; e) 

/ = -eD u W(x, I, u\ e) (x, I, u) G R 2n x R m x T m (13) 

u = DrHoix, I) + eD/W(x, I, to; e) 

where < e < 1, the dot indicates differentiation with respect to s, and J 
is the standard symplectic matrix which in our case (n = 1) is simply 

and where the subsystem, for e = 0, 

x = JD x H (x,I) (15) 

is assumed to be completely integrable, and to admit a hyperbolic fixed point 
with a homoclinic loop connecting this point to itself. 
In our case we can write the equations of motion as 



f = Mo i c sw 

8p r dp r 



(16a) 



and 
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where, again, the dot indicates differentiation with respect to s. The sub- 
system (16a) is, for e = 0, precisely as in (15) and in particular it admits 
an hyperbolic fixed point, given in (9). The subsystem (166), here writ- 
ten in terms of the spherical variables <ft, 8, should be transformed as in the 
second and third lines of (13) by introducing two action variables I (e.g., 
A = P<j>i I2 — L, the total angular momentum) with their conjugated cyclic 
variables u)\^u)i\ we shall see that it is not necessary to explicitly perform 
this transformation. 



3 The Melnikov Method 

According to |J, the general expression of the Melnikov functions for a system 
of type III reduces in our case to 

/+00 _ 
D^W{q{{s)^)ds (2 = 1,2) (17) 

-00 

where t 

q 7 (s) = (x T (s)J,J DjHoix^OJH + Uo) (18) 

is a generic homoclinic orbit of the unperturbed problem. / has to be cho- 
sen so that it defines a KAM torus (see || for more details). Actually, 
the convergence of (17) is a delicate matter: indeed the integrals converge 
only "conditionally", i.e. when the limits of integration are allowed to ap- 
proach +00 and —00 along suitable sequences T? and — TJ respectively. 
Such sequences must be chosen as in 0j, i.e. for every e sufficiently small 
we have to consider monotonely increasing sequences of real numbers, with 
j = 1, 2, . . . and hm^^T/'" = 00, such that lim^^ |g e s (T/) - q^{-Tf)\ = 
and lim^oolD^T/^O)! = lirn^ \D U {&{-T?), 0)| = 0, where q s e > u are 
trajectories of the perturbed system (e 7^ 0) in the stable and unstable man- 
ifolds (see 0). 

Moreover we recall that the existence of simultaneous zeroes of the Mel- 
nikov functions in (17) is sufficient to prove the occurrence of chaos. In fact 
it is known that the presence of simultaneous zeroes (and the periodicity of 
the perturbation) implies an infinite sequence of transversal intersections of 
the stable and unstable manifold leading to a chaotic dynamics (see e.g. [[J). 
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With 1\ = I2 = L, and dropping the superscript over I to simplify 
the notation, we then find from (17) the two conditions 



/+00 
{P4,,W} ds= 0. (19) 
-00 



and 

M L = I {L, W} ds = (20) 



where both integrals, according to (17-18), are to be evaluated along a generic 
homoclinic orbit q T Q (18). 

Apart from a rotation, which transforms the given perturbation W into 
a new W (see below for more details), we can always choose the homoclinic 
orbit in the plane 9 = ir/2, therefore it is easily seen that (19-20) become 
simply 

f +oc dW 

M p ,= —da = (21) 



and _ 

M L = J P(j) — ds = (22) 

Condition (22) turns out to be the same as (21), and then we are left with 
only one condition. Hence, to find transversal intersections of the stable and 
unstable manifolds to the periodic orbit, and therefore a chaotic behavior, 
we have to find simple zeroes of the Melnikov function (21). Let us remark 
incidentally that the transversality of these intersections is not strictly nec- 
essary, indeed - to have chaos - it would be sufficient that the crossing is 
"topological" , i.e. that there is really a crossing from one side to the other 
@- 

It is interesting to remark that we have found that, in our case, only 
one Melnikov condition has to be studied. The same result would be clearly 
obtained for the same problem restricted to the plane. This shows that to 
study the occurrence of chaos in this kind of systems (i.e. with conserved 
Hamiltonian and angular momentum) it is sufficient to consider the planar 
problem. 

As already stated, to transform the conditions (19-20) into (21-22), a 
rotation is necessary, and in particular one may choose the rotation in such 
a way that the generic homoclinic orbit is transformed precisely into the 
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homoclinic orbit in the plane z = with axis coinciding with the x axis, 
as chosen in Sect. 2 (see (10-11)). This rotation is defined by the following 
Euler angles (with the conventions and notations as in [1] 



-Q, -i, -uj (23) 

where (with the language of celestial mechanics) i the inclination of the plane 
of the orbit, uj the angle of the perihelion with the line of nodes in the orbital 
plane, and Q is the longitude of the ascending node. Notice that, in terms 
of our previous variables, one has i = arccos (p^/L) and 

Q — ojqi, uj = uj 02 (24) 

which then play the role of the arbitrary "phases" ujq in eq. (18). 

Denoting by A the matrix of this rotation, the given perturbation W(x, p x ) 
(x = (x, y, z), etc.) will assume a new expression W obtained replacing x with 
Ax and p x with Ap x . It is then easy to verify that the Melnikov condition 
(21) becomes 

M(u,tt)= (25) 
/ R(s)(-C 1 (R(s),^(s),oj,Q) sin$(s) + C 2 (. . .) cos $(s)) ds = 

J — oo 



with 
Ci 



— — ) (cos uj cos Vt — cos i sin uj sin Q) + (26) 

OXi/0 



f 9W \ , o • • on ( dW \ ■ 

+ — — (— cosu;smi2 — cos « sin u; cos i 2) + I— — si 

V ox? / o V dx 3 J o 



sin % sin uj 



fdW\ 

Ci = (— — ) (sincjcosfi + cosicosa;sinr2) + 

V ox\ / o 

fdW\ . . . . ON fdW\ . . 

+ — — (— sma;sini2 + cos z cos a; cos sZ J — — — smzcosa; 

V ox 2 / V ax 3 / 

and where (dW / dxi)o means that in the derivative of the given W with 
respect to Xj one has to replace x with Ax. and finally put z = (or 6 = tt/2). 

If, e.g., i = 0, i.e. if the problem is completely planar, including the per- 
turbation, or if the perturbation is "generic", i.e. has no "preferred" direction 
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in the space (as often happens, see next Section for explicit examples), and 
therefore it is not restrictive to assume i — 0, then one gets 

^ ( d W\ ± (dW\ . ± „ (dW\ . ± (dW\ 
Ci = - — cos 0o +1-5 — sin 0o, C 2 = [-= — sin 0o -1-5—) cos0 o 
V oxi / o V 0x2 ' V ox\ / V 0x2 ' 

(27) 

where 0o = — (u + fi), the needed rotation is simply a rotation of angle 0o 
around the z— axis and (25) becomes 

r + °° dW 

M ^ = / T^T^)' + ^ n / 2 > Ri - S ^ L > °) ds = ( 28 ) 
J-00 # 

In conclusion, it is clear that verifying Melnikov conditions for the appearance 
of chaotic behaviour amounts to verifying the existence of values of u, Q for 
which (25) (or (28) in the above hypothesis) is satisfied. 



4 The Perturbations 

In this section we want to consider the motion of a relativistic charged particle 
in a gravitational field with an electromagnetic perturbation. To this end, 
we need the Hamiltonian describing a particle in an electromagnetic field, 
i.e. 

H = ^-9 a \p a - eA a ) (p b - eA b ) (29) 
and hence, for weak fields: 

H = H g ab p a Ab + higher order terms (30) 

m 

Let us remark that we neglect the effect of the electromagnetic fields back on 
the metric since direct effects of the magnetic field on a charge are generally 
very large compared to indirect gravitational effects on the mass arising from 
gravity of the field energy; we recall, incidentally, that the metric of a black 
hole immersed in a uniform magnetic field was found as an exact solution of 
the Einstein-Maxwell equations [18, 10 1. 



Let us first consider a uniform electric field in a generic direction 1; then 
Ai = (— 0, A) = (— £ h ■ 1, 0) = {—£{l\ cos sin 6 + I2 sin sin 6 + l 3 cos6*), 0), 
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where h is the unit vector in the r direction. The perturbed Hamiltonian 
(30) is then 

H = H Q - eEf-\k cos sin 9 + l 2 sin sin 9 + l 3 cos 9) (31) 

where e = e^/m<l . 

More interesting, and a little more complicated, is the case where the 
perturbation is given by a mag netic field. Let A, = (0, A) = (0, |B x r) 
where: 

A = — (k 2 z - k 3 y, k 3 x - kxZ, k x y - k 2 x), (32) 

with B = B k = B {k\, k 2 , k 3 ), where k is a unit vector in the direction of 
the magnetic field. Rewriting the vector potential in spherical coordinates 
we obtain 

A r = 

= ^2~(k 3 sin 2 9 — sin 9 cos9(ki cos + k 2 sin 0)). (33) 

Ae = ^Tj-{k 2 cos — k% sin 0) 

In both cases the direction of the perturbing field is generic, therefore no 
rotation is required (it would simply change the directions 1, k, which are 
not fixed; see the remark at the end of previous section). 
The perturbed Hamiltonian (29) can be written as 

B cr 2 

H = H \g H P<b{h sin 2 9 - sin 9 cos 9(k x cos + k 2 sin 0))1 - 

2m 

Bcv^ 

— 2m~ [(9 e0 Pe(h cos (f) - h sin (f)))] + 

B 2 e 2 r i 2 

H \g H (h sin 2 9 - sin 9 cos 9(h cos + k 2 sin 0))1 + 

8m 

+ ^^[(/ e (A; 2 cos0-A; 1 sin0))] 2 (34) 
8m 

where this time we retain also the quadratic terms in the magnetic field, in 
view of the discussion in the next section. Since we need in ([H]) (or (f2~8|)) 
quantities evaluated along the homoclinic orbit in the plane 9 = tt/2, we can 
rewrite, recalling that = r~ 2 (sin9)~ 2 and g ee = r~ 2 , the Hamiltonian 
(34) as: 

9 

H = H - ek 3 P<f) + e 2 — - [{k\ - k\) sin 20 - k x k 2 cos 20] (35) 
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where e = § « L 



5 The Melnikov Conditions 

The Melnikov integral for the perturbation produced by the uniform electric 
field can be found from (ET) (or (|2"8"|)): 

/+oo 
Ef-\h sin($(s) + O ) - Z 2 cos($(s) + O )) <te (36) 
-oo 

since 6 = ir/2. Now, using the fact that R and $ are respectively even and 
odd functions of s, we can write the integral as: 

/+oo 
Ef" 1 cos$(s) ds. (37) 
-oo 

or, defining two constants L and a such that l\ = Lcosa and I2 = Lsina, 
the Melnikov condition is 

= J X L sin(0 o -a) = (38) 

where 

/+00 
Ef^cos^is) ds (39) 
-00 

From numerical evaluations and general arguments for this type of (condi- 
tionally convergent) integrals (see also 0), we can assume that J\ 7^ 0, or 
that it vanishes for at most some isolated values of the parameters (p^, m, M) 
involved. Therefore, the Melnikov function (j37D has simple zeroes and hence, 
thanks to the periodicity of the perturbating term in the integral, there is 
a infinite sequence of transversal intersections of the asymptotic stable and 
unstable manifolds, leading as well known to a chaotic dynamics. 

In the case of the magnetic field, it is clear from (35) and (28) that the 
first-order term in e of the perturbation gives no contribution to the Melnikov 
integral: actually, it is easy to see that the perturbed Hamiltonian truncated 
at the first-order is integrable (indeed, it admits L 2 as an additional constant 
of motion). Then, if we consider the first non- vanishing contribution to the 
Melnikov integral, this gives the following condition 

/+00 
R 2 (s)cos2$(s)ds = (40) 
-00 
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which can also be written, with obvious notations, similar to the above (37- 
39), 

J 2 K sin(2(j) - 5) =0 (41) 

As before, we can conclude that chaos occurs also in this case, but - unlike 
the case of electric field - this is now a "second-order effect". On the other 
hand, the motion of test particles moving around a black hole immersed in a 
magnetic field was studied in QT0j| , where the authors, performing a numerical 



study of the orbits of a particle in Ernst space-time, presented strong evidence 
of the occurrence of chaos and nonintegrability. As the authors acknowledge, 
the numerical methods they use cannot give a rigorous proof of the noninte- 
grability or of the occurrence of chaos. Instead, the Melnikov method is able 
to detect the existence of a infinite sequence of transversal intersections of 
the stable and unstable manifolds and hence, via the Smale-Birkoff theorem, 
to manifest the equivalence to a symbolic dynamics expressed by the Smale 
horseshoe. Therefore, the Melnikov technique provides an analytic proof of 
the occurrence of chaos in the problem discussed here and in [[1(|. In par- 



ticular our analysis shows that the chaotic dynamics appears even when the 
reaction of the magnetic field on the black hole is neglected, but does not 
when the terms quadratic in the magnetic field (|29|) are neglected. 

Finally, it can be observed that, in both cases (i.e. for both the electric 
and magnetic field) whereas the component of the field on the plane of motion 
leads to a chaotic dynamics, the component normal to the plane does not. 
This behavior can be explained observing that the component normal to the 
plane of motion is a constant on such plane. Therefore it doesn't break the 
symmetry of the system on the plane of motion, and hence it doesn't lead 
to the appearance of chaos. Since the problem is spherically symmetric, the 
reasoning used for the 9 = n/2 plane can be applied to every plane for the 
origin of coordinates. Thus, given an electric or magnetic field, on each plane 
for the origin (except at most the one normal to the field) , chaos appears for 
a suitable choice of initial conditions. 
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